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Abstract 

A complete set of filters Fn for the optimal-depth n-input sorting net¬ 
work problem is such that if there exists an n-input sorting network of 
depth d then there exists one of the form C (B C for some C £ Fn. Pre¬ 
vious work on the topic presents a method for finding complete set of 
filters Rn,i and Rn,2 that consists only of networks of depths one and two 
respectively, whose outputs are minimal and representative up to permu¬ 
tation and reflection. Our main contribution is a practical approach for 
finding a complete set of filters Rn,3 containing only networks of depth 
three whose outputs are minimal and representative up to permutation 
and reflection. In previous work, we have developed a highly efficient al¬ 
gorithm for finding extremal sets ( i.e. outputs of comparator networks; 
itemsets; ) up to permutation. In this paper we present a modification 
to this algorithm that identifies the representative itemsets up to per¬ 
mutation and reflection. Hence, the presented practical approach is the 
successful combination of known theory and practice that we apply to 
the domain of sorting networks. For n < 17 , we empirically compute the 
complete set of filters i?n,2, Rn,3, Rn,2 \ w and i?™ 3 of the representative 
minimal up to permutation and reflection n-input networks, where all but 
Rn,2 are novel to this work. 

1 Introduction 

A sorting network is an abstract mathematical model designed to sort numbers 
in a predetermined sequence of comparators. A sorting network consists of n 
wires and comparators between pairs of wires such that any input of n numbers 
is sorted by the network, where one wire corresponds to one number. The two 
most common measures of sorting networks are the total number of comparators 
— Bose-Nelson’s sorting problem m — and the number of network levels, also 
referred to as depth. In this paper, we extend the knowledge on the problem 
of empirically searching for sorting networks of minimal depth by significantly 
reducing the number of candidate network levels of depth three that need to be 
considered by any algorithm. For example for n = 11, contains about 10 
thousand comparator networks in comparison to existing methods that would 
consider roughly 1.7 million eleven-input comparator networks of depth three, 
that is about 169 times fewer networks to consider in comparison to applying 
all levels. 
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2 Related Work and Contributions 


Knuth [2] showed the optimal depth sorting networks for all n < 8. He also 
presents the zero-one principle of sorting networks which states that if a com¬ 
parator network sorts all 2" binary strings of length n then it is a sorting 
network. 

Parberry [3] presented a computer assisted proof for the minimal depth of 
a nine and ten-input sorting networks. He significantly reduced network level 
candidates for the first two levels, in comparison to the naive approach, by 
exploiting symmetries of the networks (referred to as first and second normal 
form 0). 

Bundala [4] presented a computer assisted proof for the optimal depths of 
networks with eleven to sixteen (inclusive) inputs. He also managed to signif¬ 
icantly reduce the number of candidates for the second layer in comparison to 
Parberry’s approach, by considering only networks whose outputs are minimal 
representative up to permutation and reflection. Similar work for the second 
level is also presented by Michael Codish in . Bundala’s algorithm for finding 
sorting networks of optimal depth is based on a SAT encoding of the optimal 
depth sorting networks problem, which uses the set of candidate two-layer net¬ 
works as a fixed entry point. Some extra pruning techniques are presented and 
they use a state of the art “off the shelf” SAT solver to find the optimal depth 
sorting networks for all n < 16. The work presented in this paper can be used 
to speedup Bundala’s algorithm substantially, as described in our experiments 
section [5] 

We [B] presented a highly efficient practical algorithm for finding the mini¬ 
mal representative itemsets over a domain H up to a permutation of D. This 
algorithm can be applied to reduce the number of candidates for the second 
layer as described in Bundala and Codish, although [4] and [5] present an extra 
pruning method using reflection. 

2.1 Problem Statement 

The problem addressed in this paper is that of reducing the candidate networks 
of depths two and three that need to be considered when searching for sorting 
networks of optimal depth. 

2.2 Motivation 

The importance of finding the Rn ,3 and Rn ,3 \ w set of filters is easily seen 
from Bundala’s algorithm for finding sorting networks of optimal depth and is 
explicitly stated in the future work section of [1]. Bundala’s algorithm can be 
easily adapted to use prefixes of exactly three layers as an entry point to the SAT 
encoding of the problem, given that the presented algorithm [3] uses exactly two 
layers as the entry point. Hence, such a reduction of the search space would 
result in a faster such SAT-solver-based algorithm for finding sorting networks 
of optimal depth. 

It is important to mention that the techniques described by Bundala [3] and 
Codish for reducing the number of candidate networks of depth two cannot 
be easily extended for networks of depth three because they are based on finite 
case studies/identification around the second level. In other words, they provide 
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a regular expression for generating the networks of depth two whose outputs are 
minimal and representative up to permutation and reflection. Identifying all of 
the cases to derive a regular expression for the third level is a possible solution 
to the problem but could prove an immensity hard one. 

Moreover, using our approach we manage to further reduce the required set 
of two-layered networks that needs to be considered when finding optimal depth 
sorting networks. Since Bundala showed that is a sufficient to find a subset of all 
of the 2" possible inputs for which no n-input sorting network of depth d exists 
that sorts all of them when proving depth optimality. Intuitively, the more we 
reduce the initial input set, the lower the number of minimal representative up 
to permutation and reflection networks of depth two is enough to be considered. 
We present experiments showing the achieved reduction for all n-input two 
layered networks for n < 17. 

2.3 Contributions 

• Modified algorithm for finding minimal itemsets up to permutation and 
refleetion — we took an existing algorithm for finding minimal itemsets 
up to permutation which a dataset (a collection of itemsets). We present 
a modihcation which is linear (in terms of the number of itemsets) in 
time and space to find the ones which are minimal up to permutation and 
reflection. 

• Empirically find Rn ,2 \ w for all n < 17 — this is a direct improvement of 
Bundala’s technique of not considering all of the 2" inputs to determine 
that no n-input sorting network of depth d exists. We take one step 
further to hnd the minimal up to permutation and reflection itemsets 
after applying the input set reduction (described in Experiment 3 [J]). 

• Empirically find Rn ,3 and R^ ^ for all n < 14 — we experimentally evalu¬ 
ated the modified algorithm to find the three layered n-input comparator 
networks whose outputs (itemsets) are minimal up to permutation and 
reflection. The set Rn.s is generated by applying all network levels to 
i?„ 2 and then hnding the minimal up to permutation and reflection ones, 
whereas R^^ is derived by applying all levels to all itemsets in the set 
Rn .2 r w and then reducing. 

3 Background on Sorting Networks 

3.1 Formal Definition of Comparator and Sorting Net¬ 
works 

Definition 3.1. A generalized comparator is an ordered pair {i,j) sueh that 
1 i j if n. A generalized eomparator is a comparator or min-max com¬ 
parator if i < j. The values i and j are referred to as channels. A gen¬ 
eralized level L is a set of generalized comparators such that each channel is 
involved in at most one generalized comparator, formally if {a,b), {c, d) € L 
then |{a, 6, c, fi}|= 4. A generalized level is a level or min-max level if it con¬ 
sists only of (min-max) comparators. The set of all (min-max) levels is denoted 
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as Gn, as described by Bundala A generalized n-input comparator net¬ 
work is a vector {Li,L 2 , ■ ■ ■, Ld, n), where Li, L 2 , ■ ■ ■, Ld are generalized levels, 
and n is a positive integer. A generalized n-input comparator network is called 
an n-input comparator network if it consists only of (min-max) levels. Let 
C = (Li, L 2 , ■. ■, Ld, n) he an n-input comparator network, we define the size of 
C as the number of levels, i.e. \C\ = d. 

So far we have formally defined the structure of a (generalized) comparator 
network. We need to define the output of applying a comparator network to 
an input, where an input is an n-bit binary string [5]. Applying a network to 
an input permutes the input vector. Hence, for any fixed input we can define a 
permutation that models the network behaviour when applied to that particular 
input. 

Notation 3.2. Denote the set of all permutations of n elements as n„ = {tt : 
{1,2,..., n} I— > {1, 2,..., n} | tt is bijective }. Let v = { 01 , 02 , .■.) be a vector. 
Denote by Vi the i-th coordinate ofv, namely Vi = Oi. 

Definition 3.3. An input is a vector x £ {0,1}" as per Knuth’s zero-one 
principle. Denote by /„ the set of all inputs. The evaluation of a generalized 
n-input comparator network C = {Li,..., Ld, n) in channel i at level k on input 
X is the two dimensional vector ex{i,k) where: 


ex{i,k) 


' {xi,i) 
ex{i,k- 1 ) 

ex{j,k- 1 ) 

ex{i,k- 1 ) 

ex{j,k- 1 ) 

^ex{i,k- 1 ) 


k = 0 

{i,j) G Lk and ex{i,k- l)i <= ex{j,k - l)i 
{i,j) G Lk and ex{i,k- l)i > ex{j,k - l)i 
{j,i) G Lk and ex{i,k - l)i >= ex{j,k- l)i 
(j,*) G Lk and ex{i,k - l)i < ex{j,k - l)i 
otherwise 


The output of applying C to x is Vc{x) = {ex{^, d)i,..., ex{n, d)i) £ The 
permutation of the coordinates when applying C to x is Pc{x) = {ex{^, d) 2 ,.. ■, ex{n, d) 2 ) £ 


Intuitively, we say that a vector in /„ is sorted if its values are non-decreasing 
left-to-right, and a sorting network is one which sorts all possible 2" input 
vectors. More formally: 

Definition 3.4. The vector {x\,X 2 ,.. ■ ,Xn) £ {0,1}” is sorted iff Xi <= Xj+i 
for all 1 < i < n. A generalized sorting network is a generalized n-input 
comparator network for which there exists a permutation tt £ n„ such that 
7r(Vc'{a;)) is sorted for all inputs x G In- A sorting network is an n-input 
comparator network such that Vc{x) is sorted for all inputs x G In- 

Theorem 3.5. For every generalized sorting network there is a sorting network 
with the same size and depth. If the former has only min-max comparators in 
the first k levels, then the latter is identical in the first k levels. 

Proof. See Knuth [2]. □ 
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3.2 Known Overlapping Theory 

The following definitions, lemmas, theorems and proofs presented in this sec¬ 
tion are developed independently of the work by Michael Codish 13 0 and 
Daniel Bundala [4] . But since the results described in this section have already 
been published, we label them as known properties of sorting networks. 

Definition 3.6. Let the output set of a comparator networkC he Sc = {Vc'(a;)|a; G 
/„}. Let the set of all already sorted inputs Tn = {(xi,X 2 , ■ ■ ■ ,Xn) \ Xi^j = 

0, Xiy=j = 1 for 1 < j <n-|-l}. 

Definition 3.7. Let A and B be n-input comparator networks, where A = 

{Ai, A 2 , ■ ■ ■, Ad,n), B = {Bi, B 2 , ■ ■ ■, Bk,ri), and let L be a level. Define the 
concatenations A® L = {Ai ,..., Ad, L, n) and A®B = A® Bi® B 2 ® ■■■ ® B}:. 

Note that © is associative. 

Theorem 3.8. Let A, B and C be n-input comparator networks. Suppose 
that Sa C Sb and B ® C is an n-input sorting network. Then there exists a 
comparator network C with the same depth as C such that A®C' is an n-input 
sorting network. 

Proof. See proof of the more general Theorem 13.101 □ 

Knuth [5] has shown that comparator networks are just as powerful as gener¬ 
alized comparator networks. He shows that the group of generalized comparator 
networks is closed under permutation. Intuitively, we would like to strengthen 
the result of Theorem 13.81 by considering permutations of output sets. Before 
we present this result, we need the following lemma to prove it. 

Lemma 3.9. Let tt G n„, x G In, and C be a comparator network such that 
Vc'(7r(a;)) is sorted. Then 7r(I4—i(c)(3;)) sorted, where Tr~^{C) is a generalized 
comparator network. 

Proof. Let x = {xi,X 2 , ...,Xn) and Pc{tt{x)) = {pi,P 2 , ■ ■ ■,Pn). Applying 

to the equality yields that P„-i(C)(7r"i(7r(x))) = P^-^c){x) = (7r-i(p^-i(i)),..., 7r-i(p^-i(„))). 

Applying tt to the equality yields Tr{P^-i(^c}ix)) = (pi,P2, ■ ■ ■ ,Pn) = Pc{t^{x)). 

From Definition l3.3l of the functions Vc{x) and Pc{x), we now have that 7r(T4-i((;7)(x)) = 
VciTrix)). From the hypothesis we know that Vc(7r(a;)) is sorted, hence we con¬ 
clude that 7r(t4-i((7)(a;)) = Vc{'x{x)) is sorted. □ 

Theorem 13.81 tells us that if we can extend the comparator network B to 
a sorting network by appending I levels to it then we can extend any network 
A such that Sa Sb hy appending I levels to it. We now extend this result 
by weakening the constraint Sa C Sb- We show that it is enough to find one 
permutation tt G n„ such that tt{Sa) C Sb to claim that if we can extend the 
comparator network P to a sorting network by appending I levels to it then we 
can extend any network A by appending I levels to it. 

Theorem 3.10. Let A, B and C he n-input comparator networks, and tt G n„ 
such that tt^Sa) Q Sb and B ® C is an n-input sorting network. Then there 
exists a comparator network C with the same depth as C such that A® C is 
an n-input sorting network. 
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Proof. From the hypothesis we know that there exists C such that B 0 C is an 
n-input sorting network. From tt^Sa) Q Sb we deduce that Vc'(7r(a:)) is sorted 
for all X G Sa because Tr{x) € Sb- Applying Lemmato all x G Sa and C we 
deduce that 7r(T4-i(C)(^)) is sorted. Hence H©7r“^(C') is a generalized n-input 
sorting network of depth k. Finally we apply Theorem 13.51 to the generalized 
sorting network H©7r“^(C') to show that there exists a comparator network C 
with the same depth as C such that H © C" is an n-input sorting network. □ 

Definition 3.11. Let X be a set of output sets of n-input comparator networks. 
Define the set of all minimal representative output sets up to permutation of X 
as MinPi{X) = {Sa \ Sa G X : $ Sb G X,Tr G H^, : B < AjTt^Sb) C S'a}, 
where by B < A we denote the lexicographic order of networks, as described 
by Barberry Let the set of all output sets of n-input comparator networks 
of depth d be defined as Gn,d- Let the set of all minimal representative output 
sets of n-input comparator networks of depth d up to permutation be defined as 
Sn,d = MinPl{Gn,d). 

Definition 3.12. The set Xn of n-input comparator networks is a complete 
set of filters iff for any n-input sorting network of depth d there exists one of 
the form G : C' of depth d for some C G A„. We would also denote the set 
of all complete sets of filters of that contain only n-input comparator networks 
with exactly i levels as Fni = {A„ | is a complete set of filters and C G 
Xn ^ Id = f}. 

3.3 Known Non-Overlapping Theory 

Definition 3.13. Let x = (xi,X 2 , • ■. ,x„) S then x^ = {xf[, Xn-i, • ■ ■ ,xa) 
where Xi G {0,1} and 0 = 1 and 1 = 0. Let L be a level the its reflection 
L^ = {(n — j +1, n — i +1) | {i,j)}. Let G = {Li, L 2 ,..., Ld, n) be a comparator 
network then its reflection G^ = {Lf^, L 2 , ■ • ■, L^, n). 

Lemma 3.14. Let G be a comparator network then x G Sc G Squ. 

Proof. Refer to the proof of Lemma 8 in ID by Michael Codish. □ 

Lemma 3.15. Let Rn,i be the set of minimal representative up to permutation 
and reflection itemsets within Gn. Then Rn,i G Fny- 

Proof. Refer to section 4.2 in a. □ 

Definition 3.16. The set of inputs B \ w = {{xi,... ,Xn) G In \ Xi = ... Xi = 
0, Xn-r+l = ■ ■ ■ = Xn = l,l + r = w}. 

Remark 3.17. Bundala noted that if there exists an input set B G In for 
which no n-input sorting network of depth d exists then there does not exist an 
n-input sorting network. The input sets considered by his method are of the 
form B \ w, 

Note, that w is strongly dependant on n. 
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4 Algorithm Modification 

In this section we present a modification to the algorithm for finding minimal 
itemsets up to permutation that allows us to find the minimal up to permutation 
and reflection as per Definition 13.131 and Lemma 13.151 The pseudo code of the 
modified version is presented in Algorithm [TJ 

4.1 Detailed Description 

The first thing the algorithm does is given a dataset F it makes sure that for 
every itemset Fi G F its reflection Ff^ is also in F, if this is not the case 
then we add it to the input dataset F. We assume, that the version of the 
algorithm for finding minimal sets up to permutation returns an array of 1^1 
integers subset^of such that if the itemset Fi is minimal up to permutation 
we have subset_of[i] = i; otherwise there exists a permutation tt G n„ such 
that TT{F[subset-of[i]]) C A[z]. This array gives us detailed information about 
which itemset is a subset (up to permutation) of another. We use this extra 
information to remove the ones which are non-minimal up to reflection from the 
ones which are minimal only up to permutation. We do this by iterating through 
the list of itemsets that are minimal up to permutation and for each itemset at 
index i we find the index reflect[i] of its reflected itemset within F. Then we 
traverse the subset^of[reflect[i]] until we reach the index of an itemset that is 
minimal up to permutation. If the index at which we arrived is smaller than 
the index i then we mark the minimal up to permutation itemset Fi as non- 
minimal up to reflection because we choose the lexicographically (index-wise) 
smallest itemset as representative up to reflection. This modified algorithm is 
easily proven to return the minimal up to permutation and reflection itemsets 
within F by using Definition 13.131 and Lemma 13.151 

4.2 Complexity Analysis 

The worst case time and space complexity of this modified algorithm are the 
same as the unmodified version because all we do is add an extra 0(|A|) time 
and space to the existing approach. Hence the worst time and space complexity 
are 0( ’'^”'pll~^ll ) (using P parallel threads) and 0(||F|| +r X n^) respectively. 

5 Experimental Evaluation 

We have summarized the results of our experiments in Figure [TJ We show the 
sizes of the sets Rn.i, Rn,2 for all n < 17 and Rn,3 for all n < 14. The sizes 
of Rn.i and Rn,2 match exactly to the ones presented by Bundala [4] whereas 
the set Rn,3 is novel to the work presented in this paper. The row 
presents the expected speedup of Bundala’s algorithm if the first three layers 
are to be fixed rather than only the first two — as is described in [Tj. This is 
technique is most useful when using Bundala’s method to find the satisfiable 
instances when checking if an n-input sorting network of depth d -|- 1 exists; i.e. 
for n = 12 we do expect Bundala’s program to execute about 60 times faster to 
determine that there exists a sorting network of depth d -f 1 = 8. 
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n 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

d 

4 

4 

5 

5 

6 

6 

7 

7 

8 

8 

8 

8 

u 

2 

2 

2 

3 

3 

4 

4 

5 

3 

4 

7 

7 

|G„l 

26 

76 

232 

764 

2 620 

9496 

35 696 

140 152 

568 504 

2 390480 

10349536 

46 206 736 

IRn.ll 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

|Rn,2| 

4 

5 

8 

12 

22 

21 

48 

50 

117 

94 

262 

211 

|R,n 2 t ^1 

3 

2 

3 

6 

13 

12 

20 

24 

103 

66 

83 

200 

1 |Rr,.2l 1 
L|Rnaftj|J 

1.33 

2.5 

2.67 

2.00 

1.69 

1.75 

2.4 

2.08 

1.14 

1.42 

3.16 

1.06 

|Rn,3| 

4 

4 

52 

38 

1554 

3169 

55 722 

117517 





IRt.al 

4 

4 

27 

55 

685 

971 

12 025 

38758 

2 403835 




1 |Rn,2|«|G„| 1 

L IRn.al J 

26 

95 

35.69 

241.26 

37.09 

62.93 

30.75 

59.63 





1 |R-n.2l«|Gn| 1 

L IR-li.al -1 

26 

95 

68.74 

166.69 

84.15 

205.37 

142.49 

180.80 

27.67 





Figure 1: Experimental evaluation summary presenting the sizes of the number 
of networks whose outputs are minimal and representative up to permutation 
and reflection of depths one, two and three for n < 17. The rows [ 

and J demonstrates the expected speedup of the existing [i] algo¬ 

rithm for finding sorting networks of optimal depth by fixing the first three 
layers, rather than only the first two. The row |Rn ,2 \ w| shows the necessary 
second layer networks that Bundala’s approach needs to consider when proving 
that no n-input sorting network of depth d exists. 


Bundala [1] noted that if there exists an set of inputs B € In for which no 
n-input sorting network of depth d exists then there does not exist an n-input 
sorting network. The input sets considered by his method are of the form B \ w, 
recall from Definition 13.161 In Figure [T] we present the row |Rn ,2 \ w| which 
presents the minimal up to permutation and reflection outputs of depth two 
that are restricted to be of that certain form. We use these when proving that 
an n-input sorting network of depth d does not exist; i.e. Bundala’s approach 
needs to consider only |Rn ,2 t <^-’1 set of fixed two-layered networks. Moreover, 
when generating the third layer is achieved by applying all levels to the 
itemsets from |Rn .2 f w|. Hence, adapting the Bundala’s algorithm to fix the 
first three layers, we would expect a speedup factor of [ J for the 

unsatisfied instances; i.e. for n = 12 we expect the modified version of his 
algorithm to find that no sorting network of depth seven exists about 180 times 
faster by fixing the first three layers with w = 5. 

5.1 Environment Setup 

In all of the conducted experiments we used a computer with four Intel Xeon 
CPU E7- 4820 processors. Each CPU has 8 cores clocked at 2.00GHz, equipped 
with SMB of third level cache and 128GB of main memory. Note that our 
experiments investigate the case when the entire data structure fits in main 
memory. 

5.2 Implementation Verification 

The correctness of our program was verified by calculating Rn ,2 and comparing 
to existing results [5] [4] (referred to as i?„). We have verified this for all n < 17. 





























6 Conclusion and Future Work 


This paper has extended the knowledge of the structure of comparator networks 
when constrained to the problem of finding minimal depth sorting networks. The 
current state of the art algorithm for finding optimal depth sorting networks 
fixes the first two layers, formulates the problem as a SAT encoding and then 
uses an existing SAT solver to find the answer. Using the work presented in 
this paper, we can fix the first three layers of a comparator network and then 
construct the SAT encoding. In the presented experiments we managed to find 
the three layer networks for all n < 13, where for n = 12 we do expect Bundala’s 
algorithm to execute around 180 times faster when the first three layers are fixed 
in comparison to when only the first two are fixed. 

For future work, we would like to improve the memory usage of the algorithm 
for hnding minimal itemsets up to permutation and reflection, as currently it 
requires the whole dataset to fit into main memory. This is the primary reason 
why we do not present results for any n > 17, as we had access to a machine 
with only 128G'i3 of main memory. 
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Algorithm 1: Pseudo code for finding the minimal up to permutation 
itemsets M of the input dataset F = {Pq, , -Fr-i} using T threads, 

where every Fi € F is an itemset over the domain D = {di,d 2 ,... ,dn}- 
We present a subroutine Find-Min-Rep-Perm which identifies the minimal 
representative itemsets of F using T parallel threads. It is important to 
note that in the Thread-Functor subroutine the variables index and isjmin 
are passed to the by reference, meaning that they are shared between 

threads. _ 

Input: Dataset F = {Fq, Fi,..., Fr-i} over the domain 

D = {di,d 2 , ■ ■ ■ ,dn} and the degree of parallelism P 
Output: The minimal itemsets within the dataset F up to permutation 
of D. i.e. MinTr{F) 

/* We add the reflections of the itemsets that are missing 
from the dataset F. */ 

1 F i — F\JF^; 

/* Remember the reflect indexes. */ 

for i i — 0 to |F| do 

2 reflect[i\ i — i' such that F[i]^ = F[i']; 

end 

3 subset-of <— Find — Min — Rep — Perm{F, P); 

4 for i i — 0 to |F| do 

5 isjmin 4 >i[i] <— false; 

/* Check if F[f] is minimal over F up to permutation. 

*/ 

6 if subset JO f[i] = i then 

/* We set the itemset Fi as minimal up to permutation. 

*/ 

7 isjminjpi[i] <— true; 

/* We start with the reflection of F[i] and work our way 
following the subsetjof path to a minimal up to 
permutation itemset. */ 

item i — r&flect[i] while C -of[item] ^ item do 
9 item <— subset-of[item]; 

end 

10 isjminjrefl[i] <— true; 

11 if item < i then 

/* We set the itemset Fi as non-minimal up to 
reflection, because we choose the 
lexicographically smallest (i.e the one with a 
least index in the dataset) to be representative 
up to reflection. */ 

12 isjminjrefl[i] i — false; 
end 

end 

end 

13 return {Fi G F \ isjmin-pi[i] & isjminjrefl[i]}; 





